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Abstract 

The Ising model in two dimensions with the special boundary conditions of Brascamp and 
Kunz is analysed. Leading and sub-dominant scaling behaviour of the Fisher zeroes are deter- 
mined exactly. The finite-size scaling, with corrections, of the specific heat is determined both 
at the critical and pseudocritical points. The shift exponents associated with scaling of the 
pseudocritical points are not the same as the inverse correlation length critical exponent. All 
corrections to scaling are analytic. 



1 Introduction 



Second order critical phenomena are signaled by the divergence of an appropriate second derivative 
of the free energy along with the correlation length. For temperature driven phase transitions, 
this divergence is in the specific heat. Of central interest in the study of such phenomena is the 
determination of the critical exponents which characterize these divergences. Here, the concept of 
universality plays a fundamental role. The universality hypothesis asserts that critical behaviour 
is determined solely by the number of space (or space time) dimensions and by the symmetry 
properties of the order parameters of the model. The universality class is thus labeled by critical 
exponents describing the singular behaviour of thermodynamic functions in the infinite volume 
limit. In finite systems the counterparts of these singularities are smooth peaks the shapes of 
which depend on the critical exponents. Finite-size scaling (FSS) is a well established technique for 
the numerical or analytical extraction of these exponents from finite volume analyses |Q, ^, ^, ^] . 

In particular, let Cl(/3) be the specific heat at inverse temperature P for a system of linear extent 
L. FSS of the specific heat is characterized by (i) the location of its peak, Pl, (ii) its height Cl{Pl) 
and (in) its value at the infinite volume critical point Cl{Pc)- The position of the specific heat 
peak, Pl, is a pseudocritical point which approaches /3c as L ^ cxd in a manner dictated by the 
shift exponent A, 

\Pl-/3c\-L~^ ■ (1.1) 

Finite-size scaling theory also gives that if the specific heat divergence is of a power-law type in 
the thermodynamic limit, Cl ~ |/3 — /3c|~") then its peak behaves with L as 

CLm-^L"^" , (1.2) 

in which u is the correlation length critical exponent. When a = 0, which is the case in the Ising 
model in two dimensions, this behaviour is modified to 

CUPl)^ In L . (1.3) 

In most models the shift exponent A coincides with l/i/, but this is not a direct conclusion of finite- 
size scaling and is not always true. If the shift exponent obeys A > 1/z/ then ( |1.2| ) also holds if the 
pseudocritical temperature (3l is replaced by the critical one (3c- If, on the other hand, A < 
one has Cl{(3c) ~ L^°' in the power-law case 0. 

The leading FSS behaviour of a wide range of models is by now well understood. Recently, at- 
tention has focused on the determination of corrections to scaling [^, ^, |^, ^, ^, |l^, |l^, [l^ . These 
corrections may arise from irrelevant scaling fields or be analytic in L~^. Typically, numerically 
based or experimental studies involve systems of limited size where these corrections to leading FSS 
behaviour cannot be dismissed. A better knowledge of universal sub-dominant behaviour would 
therefore be of great benefit in FSS extrapolation procedures jl3, |l^, |l6|, 17, [l^, 19, 20, 21, p^ . 



The Ising model in two dimensions is the simplest statistical physics model displaying critical 
behaviour. Although solved in the absence of an external field p3], it remains a very useful testing 



ground in which new techniques can be explored in the hope of eventual application to other, less 
understood, models in both statistical physics and in lattice field theory. 

Exploiting the exactly known partition function of the two dimensional Ising model on finite lattices 
with toroidal boundary conditions Ferdinand and Fisher [14| analytically determined the 



specific heat FSS to order L^^. At the infinite volume critical point this was recently extended 
to order L^"^ by Izmailian and Hu |jl^ and independently by Salas ||2^. It was found that only 
integer powers of L^^ occur, with no logarithmic modifications (except of course for the leading 
logarithmic term), i.e., 

Cl(/?c) = Coo lnL + Co + ^-4 . (1.4) 

k=i ^ 

For the toroidal lattice, the coefficients Cqo, Cq, Ci, C2 and C3 have been determined explicitly p3. 



19| , |20| , Ferdinand and Fisher also determined the behaviour of the specific heat pseudocritical 
point, finding A = 1 = (except for special values of the ratio of the lengths of the lattice edges, 
in which case pseudocritical scaling was found to be of the form InL). 

The specific heat of the Ising model has also recently been studied numerically on two dimensional 
lattices with other boundary conditions in Refs. For lattices with spherical topology, the 

correlation length and shift exponents were found to he u = 1.00 it 0.06 and A = 1.745 it 0.015, 
significantly away from 1/u [^]. This is compatible with an earlier study reporting A ~ 1.8 
Therefore the FSS of the specific heat pseudocritical point does not match the correlation length 
scaling behaviour. This is in contrast to the situation with toroidal boundary conditions, where 
A = 1/z/ = 1 |14]. On the other hand, it was established in B that the critical properties on such 



a lattice are the same as for the torus. 

The finite size behaviour of the specific heat is related to that of the complex temperature zeroes 
of the partition function, the so-called Fisher zeroes |2^. Indeed, the FSS of the latter provides 
further information on the critical exponents of the model. The leading FSS behaviour of the 
imaginary part of a Fisher zero is [p^ ] 

lmzj{L) L-'^/" , (1.5) 
while the real part of the lowest zero may be viewed as another pseudocritical point, scaling as 

|/3,-Rezi(L)| ~L-^— . (1.6) 

Hoelbling and Lang used a variety of cumulants as well as Fisher zeroes to study universality of 
the Ising model on sphere-like lattices They reported that the imaginary part of the first zero 
is an optimal quantity to determine the exponent ly which is in perfect agreement with unity for 
all lattices studied. Regarding the FSS of the pseudocritical point, the numerical results on a 



torus are in agreement with A = l/v = 1 |14|. The pseudocritical point (from specific heat, the 
real part of the lowest zero as well as from two other types of cumulant) is not, however, of the 
Ferdinand-Fisher type. I.e., the amplitude of any 0{1/L) contribution is compatible with zero. 
While a shift exponent A = 1.76(7) is consistent with their numerical results, Hoelbling and Lang 
point out that this is not stringent and pseudocritical FSS of the form L^^ log L or even are 



also compatible with their data. Thus, while a possible leading term almost or completely 
vanishes and subleading terms are dominant, the precise nature of these corrections could not be 
unambiguously decided. In any case, FSS of the position of the specific heat peak does not accord 
with the correlation length exponent for spherical lattices and the thermodynamic limit is achieved 
faster there than on a toroidal lattice. 



In another recent study |Q involving Fisher zeroes, Beale's [15| exact distribution function for 
the energy of the two dimensional Ising model was exploited to obtain the exact zeroes for square 
periodic lattices up to size L = 64. The FSS analysis in yielded a value for the correlation length 
critical exponent, v, which appeared to approach the exact value (unity) as the thermodynamic 
limit is approached. Small lattices appeared to yield a correction-to-scaling exponent in broad 



agreement with early estimates (cj ~ 1.8 pqj). However, closer to the thermodynamic limit, these 



corrections appeared to be analytic with to = 1. On the other hand, in |17|, the scaling behaviour 
of the susceptibility and related quantities was considered and evidence for uj = 1.75 or, possibly, 
2 was presented. 

In the light of these recent analyses, we wish to present analytic results which may clarify the 
situation. To this end, we have selected the Ising model with special boundary conditions due 
to Brascamp and Kunz |27]. These boundary conditions permit an analytical approach to the 



determination of a number of thermodynamic quantities. Recently, Lu and Wu exploited this 



fact to determine the density of Fisher zeroes in the thermodynamic limit [28|. In this paper, we 
take a complimentary approach, exploiting FSS behaviour (i) to determine critical exponents, (ii) 
to determine corrections to leading scaling and (iii) to gain experience in the hope of eventual 
application to other, less transparent scenarios. The rest of this paper is organised as follows. In 
Section 2 the Brascamp-Kunz boundary conditions are introduced and the exact FSS of the Fisher 
zeroes calculated. The specific heat and its pseudocritical point are analysed in Section 3. Our 
conclusions are contained in Section 4 and the Appendix contains some calculations of relevance 
to the specific heat analysis of Section 3. 



2 The Fisher Zeroes for Brascamp-Kunz Boundary Conditions 

Brascamp and Kunz introduced special boundary conditions, for which the Fisher zeroes are known 
for any finite size lattice ||2^ . They considered a regular lattice with M sites in the x direction and 
2N sites in the y direction. The special boundary conditions are periodic in the x direction and 

Ising spins fixed to ... + + + ... and ...H 1 1 — ... along the edges in the y direction. For 

such a lattice, the Ising partition function can be rewritten as 

N M 

Zm,2N = 2'*^^ n n [l + - ^(cos 0^ + cos 0,)] , (2.1) 

j=lj=l 

where z = sinh2/3, Oi = (2i — l)7r/2A^ and (f>j = jiT/{M + 1) and where /3 = l/ksT is the inverse 



temperature. The multiplicative form of ( |2.lD is of central importance to this paper. 

Brascamp and Kunz showed that the zeroes of the partition function (^) are located on the unit 



circle in the complex z plane (so that the critical point is z = = 1). These are 



where 



= cos 



Zij = exp {ioLij) , 
^ / cos 9i + cos 4>j 



(2.2) 



(2.3) 



Setting 2N = aM and using a computer algebra system such as Maple, one may expand (|2.3| ) in 
M to determine FSS of any zero to any desired order. Indeed, the first few terms in the expansion 
for the first zero which is the one of primary interest, are 



vro" 



-M- 



V2a 
7rV2 



V2VI + 0-2 



4 7r\/2 \ a" , , _ 9 „ /l^ vr 



(1 + ^2)5/2 ]_4 

-M 

+0 



{l + a^)i3 + 2a^ 



vr 



4\2 



960-3 



(1 - a^) 



ail + .2)5/. I + + + + - + 



(2.4) 



Separating out the real and imaginary parts of the FSS of the first zero yields 



Rezii = cos ail = 



+ M-'— + M~*— I -3 + — f 1 + A 



12 



+ 0[M 



(2.5) 



and 



Imzii = sin ail 



+M- 



TT 



^/2 



2\2 



^(1 + ^2)5/2 



-7r2(l + a'^f (5 + + ba") + 24a^(l + a^){3, + 2^^ 

[-^ (4 + + 2^4) + ^(1 + ^')(1 + 3^2) (7 + 5^2)1 \+0 
4 ^ ^96 J 



(2.6) 



respectively. From the leading term in ( p.6| ) and from ( |1.5| ) one has, indeed, that the correlation 
length critical exponent is 1. Note, however, from ( p.5[ ) that the leading FSS behaviour of the 
pseudocritical point in the form of the real part of the lowest zero is 



Rez 



11 



1 - Rezii ~ M 



(2.7) 



giving shift exponent Azero = 2. This value is consistent with the numerical results of ||^] for spherical 
lattices. One further notes that all corrections are powers of and in this sense entirely analytic. 



Setting fj = 1 gives the FSS behaviour of the first zero for a square M x M lattice to be 



an = vrM"^ H M'^ 

2 8 16 



(2i 



which, for the real and imaginary parts separately is 
Rezii = cos ail = 1 
Imzii = sinaii = vrM^"*^ 



2 2 
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(2.9) 



vr 

+ -|vr- 



11 



M"^ + O ( M 



(2.10) 



In summary, we have observed that for the first zero, the shift exponent is not 1/z^ and that all 
corrections are analytic. Expansion of higher zeroes yields the same result. 

3 The Specific Heat 



In terms of the variable z = sinh 2(3, the specific heat is 



C = 



[l + z' 



52lnZ d\nZ' 



dz^ 



+ z- 



dz 



(3.1) 



where V is the volume of the system. The singular behaviour comes from the first term only, the 
second being entirely regular. Thus we may split (^]^) into singular and regular parts and retain 
only the former. From ( |2.1| ), this singular part of the specific heat for an M x 2N lattice is (up to 
the factor (1 + z2)4A;ij/32) 



smg. 
M,2N 



(^) 



1 



N M 



2MiV^^ 1 1 + z2 - zicosOi + cos(/)j) [1 + ^2 - z(cos6lj + cos(/)j)]^ 
An alternative approach is to write the partition function as 

N M 

Zm,2N = M^)YIYI{^ - Zij){z - Zij) , 
i=lj=l 



[2z — (cos 9i + cos (l)j)y 



(3.2) 



(3.3) 



where the product is taken over all zeroes. The non-vanishing function A{z) contributes only to the 



regular part of the partition function, and from (2.2) and (2.3), the product term leads precisely 



to the singular expression (|3.2|) . 



It is straightforward to perform the i- summation first in (|3.2| ). Indeed, ( [A.3| ) and ( A. 4 ) yield the 
exact result, which is conveniently expressed as 



I I M , ... A/ 



i=i 



z 



1 

2^ 



where 



9j{z) 



tanh (^N cosh ^ {1/z + z — cos (pj ) 

^{1/z + Z — COS (/>j)2 — 1 



(3.4) 



(3.5) 



and where g^p {z) = dl' gj{z) / dz^ . 



It is convenient at this stage to introduce the sums 



, M 

Su = liY.9f\l) , (3.6) 



and to make the observation that 



fff (1) = , (3.7) 

gf\l) = -3gf{l) . (3.8) 

Specific Heat at the Critical Point: At the infinite volume critical temperature, z = Zc = 1, 
the specific heat is, from (p.4[) , 

Ct,%il) = So-l . (3.9) 
The sum 5*0 is given in ( [A. 25 ) in terms of the ratio p defined through 

N = p{M + l) . (3.10) 

This gives for the critical specific heat, 

cSa)^^^^fi + i)-„.|i.-|,.-|,.o(-I,) . (3.n) 



vr 



where 



1 / 31n2 , . x\ 1 , ^ 

CO = -i^jE + ^-ln7r-2Wi{p)j-- , (3.12) 

1 / 31n2 vr , A , , 

ci = -i^jE + ^-ln7r + l-^-2Wi{p)j , (3.13) 

C2 = ^ + I^ + i^2(p)-^M^3(p) , (3.14) 

C3 = -^-T^-i^^(/') + ^^3(.) , (3.15) 

where je ^ 0.577 215 6649 ... is the Euler-Mascheroni constant and where the functions Wk{p) are 
given in the appendix. Typical values of the constants cq-cs are given in Table |l[ 



Table 1: Values of the coefficients cq-cs for various values of the ratio p = N/{M + 1). 



p 


1/2 


1 


2 


Co 


-0.376 674 2314... 


-0.350 879 7312... 


-0.349 694 204 8... 


Cl 


0.264 858 959 5... 


0.290 653 459 7... 


0.291838 986 0... 


C2 


0.125 8961381... 


0.175 784 346 0... 


0.180 950 438 8... 


C3 


-0.019 792 842 7... 


-0.069 681050 6... 


-0.074 847 143 4... 



So for the critical specific heat, apart from the In M/M term, the FSS is qualitatively the same as 



(but quantitatively different to) that of the torus topology (see ( |1.4| )and |14,[T^, ^]). We note that 
with the Brascamp-Kunz boundary conditions, it is far easier to extract the FSS behaviour. Indeed, 
determination of ©(l/M^) and higher terms proceeds is a similarly straightforward manner. 

Specific Heat near the Critical Point: The pseudocritical point, z^/2n°^ value of the 

temperature at which the specific heat has its maximum for a finite M x 2N lattice. One can 
determine this quantity as the point where the derivative of Cm,2n{z) vanishes. 

The specific heat is given in (|3.4| ). This may now be expanded about the critical point z = 1. This 
Taylor expansion simplifies using ( |3.7D and (|3.8| ). Indeed, 

^52 + 6S'o-^ 



CmM^) = -^o - ^ + (^ - 1) [-3^0 + 1] + (^ - 1)' 



+(z - 1)3 [-5^2 - 105o + 2] + O ((z - 1)^) . (3.16) 
The sums Sq and ^2 are given in ( [A. 25 ) and ( A.28| ) respectively. 



Prom ( |3.16| ), the first derivative of the specific heat on a finite lattice near the infinite volume 
critical point is 

M,2N\ { ^ ;l - 35o + 3(z - 1) [^2 + 45o - 1] + 3(z - 1)2 [-5^2-1050 + 2] + 0((z- 1)3) . (3.17) 



dz 

Noting that, while S2 is O (M^), 5*0 is O (InM), one sees that this derivative vanishes when 

35o - 1 , [5^2 + WSq - 2] [3So - if , 3X 
^-' = 3[52 + 45o-l]+ 9[52 + 45o-lf + ^ (^^ " ) ' ^''''^ 

Expansion of ( p. 181 ) now gives the FSS of the pseudocritical point to be 

pseudo ,^ InM 62 InM , 63 , f {lnMf \ , ^InM^ , .0 M ^ ,0 iq, 



^M,2iV --^"i^2 ^M2^"3^3 ^M3^^\, M4 j ^ ^ \ J ^ ^ 

where 



7r2 



2(C(3)-2iy4(/5)-47rpI^5(p)) ' ^ ■ ^ 

TT^ 67E + 91n2-61n7r-27r- 12PFi(p) 
' 12 C{^)-2W^{p)-47rpW,{p) ' ^ ■ ' 

" m-2W,{p)-AnpW5{p) ' ^^-^^^ 
, ^ ^2 1 - 27j; - 3 In 2 + 2 In TT + TT - 7r/(4V2) + 4I^i (p) 
' 2 C(3) - 21^4 (p) - ^T^pW^ip) ■ ^ ' ' 

Again, the functions Wk{p) are given in the appendix. Typical values of the coefficients are given in 
Table ^. Thus there is no leading 1 /M or In M/M term in the FSS of the specific heat pseudocritical 
point and the specific heat shift exponent does not coincide with 1/zv = 1. This is consistent with 



the result (2^) for the finite-size scaling of the real part of the Fisher zeroes. 



Table 2: Values of the coefficients 02-4 for various values of the ratio p = N/{M + 1). 



p 


1/2 


1 


2 


0,2 


-5.696 643 550 4... 


-4.200 054 895 9... 


-4.105 6215721... 


03 


11.393 287100 7... 


8.400109 7917... 


8.211243 1441... 


62 


3.758 407422 6... 


2.430 665 892 0... 


2.360 724 066 8... 


63 


-16.015 279 5172... 


-11.127129 8521... 


-10.846 367 0579... 


c[, 


-0.376 674 233 5... 


-0.350 879 733 3... 


-0.349 694 207 0... 




0.264 858 9574... 


0.290 653 45 76... 


0.291838 983 9... 




1.309 8371214... 


0.847424 969 7... 


0.829 066 715 2... 


c^, 


-6.557959 3741... 


-4.204 047 250 7... 


-4.086 049 903 8... 




-3.589 014 676 5... 


-2.321 1149405... 


-2.254 325 418 5... 




11.7044504713... 


8.304 511256 4. . . 


8.103192 065 9... 




2.719 946 882 3... 


2.005 378 462 3... 


1.960 289 872 9... 


4 


-2.719 946 882 3... 


-2.005 378 462 3... 


-1.960 289 872 9... 



Inserting ( |3.19| ) for the pseudocritical point into ( |3.16| ) gives the FSS behaviour of the peak of the 
specific heat. This is 

(lnM)2 InM .^( {\^Mf \ 



where 



c'o = Co , (3.25) 

c'l = Cl , (3.26) 

= ^ i n27^ 

2 4 C(3)-2t^4(p)-47r/,W5(/3) ' ^ ■ ^ 
^ TT 67j; + 91n2 - GluTT - 277 - 12Vri(p) 

' 4 C(3)-2t^4(p)-47rpW5(p) ' ^ ■ ^ 
, _ vr (67i;+91n2-61n7r-2^-12l^i(p))2 

"' + 48 C(3)-2H^4(p)-4^pH^5(/>) ' ^ 

4 = -4 , (3.30) 
7r2 /„ 3 



2(C(3) - m^{p) - ^T^pW^{p)) V 4^ 

3 /' 3 In 2 \ \ 

+- fl-7s-^ + lnvr + 2Tyi(/))jj , (3.31) 

f 3 / 31n2 , , ,\2 



4(C(3) - 2T^4(/o) - 47rpTy5(p)) \ vrH ' 2 
1/ 3 In 2, .x\/6 3\ 2 11 

Again, typical values for the coefficients are compiled in Table ||. 



One remarks that, up to 0(1/M), this is quantitatively the same at the critical specific heat scaling 
of ( p.ll ). One further remarks that the higher order structure of ( 3.24| ) differs to that at the critical 
point as given by ( p.ll ) in that there are logarithmic modifications to the subdominant terms. 



4 Conclusions 

We have derived exact expressions for the finite-size scaling of the Fisher zeroes for the Ising model 
with Brascamp-Kunz boundary conditions. The shift exponent, characterizing the scaling of the 
corresponding pseudocritical point is A = 2 and is not the same as the correlation length critical 
exponent u = 1. This is consistent with numerical results for lattices with spherical topology |^. 
Corrections to FSS can also be exactly determined and are found to be analytic. This is consistent 
with the large lattice numerical calculations of Q for toroidal lattices. 

A similar analysis applied to specific heat at criticality yields results qualitatively similar to those 
on a torus [|l4| , |l9| , |20| , 23|. These results complement the finite size scaling of the zeroes. I.e., apart 



from the leading logarithm, only analytic corrections appear. 

The finite size scaling of the pseudocritical point, determined from the specific heat maximum is 
governed by a shift exponent A = 2 with logarithmic corrections. This is again compatible with 
previous numerical results (5|, ^, ^. Finally, the first few terms for the FSS of the specific heat peak 
are derived. Again all corrections are analytic, but in contrast to the specific heat at the critical 
point, here they include logarithms. 

Acknowledgements: R.K. would like to thank Wolfgang Grill for his hospitality during an 
extended stay at Leipzig University in the framework of the International Physics Studies Program, 
where this work was initiated. 



A Appendix 

Consider, firstly, the sum 

1 ^ 1 

-Y r , (A.l) 

iV^(Z-cos0,)' 

where 9i = {2i — \)'k/2N , is a positive integer and Z > 1 is independent of n. We firstly treat 



the case k = 1. We follow a similar calculation in [16| and construct 



Zj — COS 

Integrate J^{z) along the rectangular contour bounded by Rez = 1/2 + e, Rez = 2A^ + 1/2 + e 
and Ivnz = HR where R is some large real constant and e is a small real number which we take 
to be positive. Because of the periodic nature of the integrand, the integrals along the left and 
right edges of the contour cancel. Further, due to the cot nz term, the integrand and hence the 
full integral vanishes in the limit of infinitely large R. Now, J-^{z) has 2A^ simple poles inside the 



contour along the real axis and a further two coming from the simple zeroes of its denominator. 
The residue theorem then gives the exact result 



1^1 1 

— V = , ^ tanh (n cosh"^ Z 



(A.3) 



This result is also implicitly contained in ||2^. Differentiation of (A.3) with respect to Z gives the 
sum (A.l) with larger values of k. In particular, we also find 



1 ^ 

AT ^ 



Z tanh cosh'^ iVsech^ (iV cosh"^ Z 



iV^(Z-cos( 



(^2 _ 1)3/2 



(A.4) 



The purpose of the remainder of this appendix is to calculate the two sums Sq and 52 where 

M 



(A.5) 



where gf\z) = df^ gj{z) / dz^ , 



tanh (A^ cosh (1/z + z — cos i 



J {1/z + Z — COS (t)jY — 1 

and where = 7rj/(M + 1). It is convenient to rewrite gj{z) as 

1 

9j{z) = 



(A.6) 



[2 + m? — cos (pjY — 1 
1 



/(2 + m? — cos (pjY — 1 
where m? = 1/ z + z — 2. We consider, firstly, the sum 

1 



|exp 2A^cosh"^ (2 + m2 -cos(^j) + l} ^ , (A.7) 



M 

E 



:i J {2 + m? — cos I 



(A. 



and we are interested in the limit m ^ 0. This may be calculated by direct application of the 
Euler-Maclaurin formula |^^. In fact, the more general summation ( |A.^ ) has been calculated in 
[16|. Writing m = r]/{M + 1), that sum is given by 



M + 



M 



2 + 



\M+1 J 



COS 0j \ — 1 



In (M + 1) ^ 1 
vr 



+ 



vr 

1 1 



3 In 2 ^ f V2rj' 
7£; - Invr H h Gq 



2(M + 1)2 



1(y2-^' 



2 



TT 



1 1 In (M + 1) 
M + I4V2 ~ (M + 1)2 4^ 



vr 



37r'^ \ vr 



r/2 / 31n2 



+ 



1 ^^,=+0 



(M + l)3 64 



InM 

{M + iy 



+ 



1 



{M + iy 



, (A.9) 



where ~ 0.577215 664 9. . . is the Euler-Mascheroni constant |3C] and where the remnant func- 
tions are 

2n 



oo 



Go(a) 
Gi(a) 



n=l 

oo 

n=l 
oo 



2n 

n 



(-1)" / 2n 



n 



C(2n + l)(^|y , 
C(2n + 1) -) 



H,{a) = ^(-l)"(2n + l) 



n=l 



2n 

n 



C(2n + 3) 



2ra 



in which C,{n) is a Riemann zeta function |30], 



c(n) = E F 



k=l 



Setting ?7 = 0, one has the first component of the sum appearing in 5*0, 

M 1 



1 



M + 



1^ 



j=i a/ (2 - cos 4) jf - 1 



1 

+- 

vr 

Next, we consider 



M 

E' 



7£; — Invr 



3 In 2 



ln(M + l) 
vr 

1 1 



+ 



vr 



+ 



M + 14^/2 (M + 1)2 144 V(M + 1 



2 + 



Af+l 



COS ( 



1 



It is natural to define an aspect ratio, p, through 

N = p{M + l) , 

and introduce 



(A.IO) 
(A.ll) 
(A.12) 

(A.13) 



(A.14) 



-1 

{exp [2iVcosh"^(2-cos(/>j)] +l} =E^i('?) • (^-l^) 



(A.16) 



X,- = y'2r/2 + vrV , (A. 17) 

Yj = exp{2pXj) + l . (A. 18) 

Dominant contributions to ( |A.15 ) come from the small j terms. The sum may thus replaced by 

oo 

~ (A.19) 



where the expansion of the summand is 



hj{T]) = + 



1 ttV X, , pYj-l TT^ 



4 ,-4 



+ 



XjYj M + 1 1 24 x^.Yi 8Yi 12 Y^^ 



1 



(M + 1) 



. (A.20) 



Taking the limit of this quantity as r] ^ gives the second component of the sum appearing in Sq, 

M 



1 



M + lfr{ ./(2-cos</.,)2-l 



{ 



exp 



2iVcosh"^(2-cos(/)j) +l} ^ 
+ ^ 



-^W2{p) + ^W,ip) 



{M + iy 



(A.21) 



where 



Wi{p) 

W2{p) 



y 

^ nie^f^"^ + 1) ' 

n=l ^ ' 



E 



n 



E 



(A.22) 

(A.23) 
(A.24) 



The functions Wfc(/9) are rapidly converging sums which may be computed numerically. For ex- 
ample, at p = 1, f« 0.001865 707 7..., W2 ~ 0.001870 9561... and PFg 0.001874495 5.... 
Finally, from ([Al5|) , (|A?^), ( |A.14D and ([A.21|) , the desired result is 



InM 1 / 31n2 , , A InM 1 

5o = + - hi? + - Invr - 2Wi{p) + — — 

1 1 / 31n2 , vr 



(A.25) 



The sum 52 is 

^2 



£ 1 *^ 



2=1 



^2 ^ M l-2|exp 27Vcosh-^(l/z + z-cos0j) + l}' 
M § Jil/z + z-cos<j)j)^-l 



z=l 



(A.26) 
(A.27) 



The sums involved here are found by taking appropriate derivatives of ([A. 91) and ( A. 201) . The result 
for S2 is 



82 = [C(3) - 2^-4 (p) - ^npW^ip)] 



+M^[C(3)-2H^4(p)-47rpt^5(p)]-^ + 0(l) , (A.28) 
in which C(3) « 1.202 056 903 2 . . . and 

00 -1 

W^ip) = ^ — -T , (A.29) 

n^(exp(2p7rn) + 1) ^ ^ 

Typical values of these rapidly converging sums are W4 0.001 864 398 1 . . . and W5 0.001 861 360 1 . . . 
at /9 = 1. 
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